Magnetotransport in nearly antiferromagnetic metals 
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We present a theory of the magnetotransport in weakly disordered metals close to an antiferromag- 
netic quantum-critical point. The anisotropic scattering from critical spin fluctuations is strongly 
influenced by weak but isotropic scattering from small amounts of disorder. This leads to a large 
regime where the resistivity obeys a scaling form p = po+Ap « po+T :i ' 2 f(T/po, (p—p c )/po, B/p ), 
where po is the residual resistivity, B the magnetic field and p — p c > measures the distance from 
the quantum-critical point on the paramagnetic side of the phase diagram. Orbital effects of the 
magnetic field are most pronounced in very clean samples for not too low temperatures, where the re- 
sistivity for increasing magnetic field crosses over from a linear temperature dependence Ap ~ T^/po 
to a resistivity linear in B and independent of T and po. At higher magnetic fields Ap saturates at a 
value proportional to T 1 or T 2 / (p — p c ) . Deviations from scaling, the interplay of orbital and spin 
contributions of the magnetic field and experimental test of the spin-fluctuation model are discussed 
in detail. 



72.10.Di,75.30.Mb, 71.27.+a,75.50.Ee 



I. INTRODUCTION 

One of the most important unresolved questions in the 
field of strongly correlated metals is the stability of the 
Fermi liquid. In one dimension, the Fermi liquid is de- 
stroyed by strong quantum fluctuations. The strange 
properties of the normal phase of high-temperature 
superconductorsa as well as the observation of an jjnex- 
pected metal-insulator transition in Si-MOSFETsa and 
other systems have cast doubts on the applicability of 
Landau's Fermi-liquid paradigm in two dimensions. 

Even in three dimensions (3d) a number of heavy- 
fermion metals near an antiferromagnetic (AFM) quan- 
tum critical point (QCP) were shown to display striking 
deviations from conventional Fermi-liquid behavioraEj. 
The main aim of this paper is to develop a detailed 
transport theory which allows to decide experimentally 
whether the non-Fermi liquid behavior near the QCP in 
weakly disordered metals can be explained by a nearly 
AFM Fermi liquid. 

It is believed that the magnetic transition in heavy 
Fermion systems results from the competition between 
the Kondo screening of magnetic moments and the AFM 
correlations induced by the RKKY interaction!]. Pres- 
sure, magnetic fields or doping can influence this inter- 
play of screening and magnetism. This allows to fine-tune 
systems directly to the QCP where the strongest non- 
Fermi liquid (NFL) effects are found; a schematic phase 
diagram is shown in Fig. Fn In particular, the resistivity 
rises as a function of temperature with exponents smaller 
than 2, the specific heat coefficient diverges or shows a 
VT cusp, and the susceptibility shows anomalous correc- 
tions of the form T a , with a < 1. We stress that both 
magnetic and non-magnetic phases are heavy-Fermion 
metals which display, for example, the characteristically 
large specific heat. The discussion of this paper is re- 
stricted to the behavior in the paramagnetic phase close 



to the QCP. 

Despite the growing amount of experimental data, a 
common agreement on the origin of the observed non- 
Fermi liquid behavior is still lacking. At present it is 
unclear, whether a single mechanism is responsible for 
the observed behavior as different exponents have been 
reported in various compounds. At least three differ- 
ent theoretical scenarios have been proposed and are dis- 
cussed in literature. 

The first scenario is based on the assumption that in a 
heavy Fermion system below a scale Tk the low energy 
excitations are (heavy) quasi particles and their collective 
excitations. In this case the QCP should be in the same 
universality class as the weak-coupling spin-density-wawe 
(SDW) transition in a Fermi liquid studied by HertzEjtZI. 
More precisely, the so-called non-Fermi liquid behavior 
near the QCP is determined by the mutual interaction 
between Landau damped spin fluctuations and inelasti- 
cally scattered quasi-particles in a nearly antiferromag- 
netic Fermi liquid. 

Alternatively, one can envision a situation where the 
Kondo effects breaks down directly aX the transition, e.g. 
due to strong magnetic fluctuationsBcd. A Fermi liquid 
description of the transition is then not possible and it 
has been speculated that Zefjai fluctuations play the dom- 
inant role in this scenaridl3l23. 

Independent of the precise nature of the QCP, suffi- 
ciently strong disorder will certainly change the nature 
of the transition. Indeed, if typical fluctuations of the 
effective Neel temperature due to disorder are large com- 
pared to the distance of the quantum critical point (Har- 
ris criterion), disorder effects are important. Perhaps 
more importantly, even at some distance form the QCP, 
in the non-magnetic phase, rare configurations of impu- 
rities can lead to a small magnetically ordered regions 
which can dominate some of the thermodynamic proper- 
ties in a finite region around the QCP (Griffiths-McCoy 



singularities)!^ 

Up to now, a comprehensive theoretical description 
exists only., for the weak-coupling spin-density wave 
transitionE3il3. As this theory is above its upper 
critical dimension, essentially all low-energy properties 
can be pdetermined analytically. It was realized only 
recentlyc3~E3 that the calculation of transprail quanti- 
ties is quite subtle and earlier predictionsEJO are not 
valid in weakly disordered materials. The complications 
arise because the scattering from spin-waves is extremely 
anisotropic and effective only in small areas on the Fermi 
surface. Therefore the transport properties strongly de- 
pend both qualitatively and quantitatively on how other 
scattering mechanisms redistribute quasi particles and 
scatter them into these small regions. 

The resistivity is often the most sensitive experimental 
probe to study the QCP, and in a number of (pressure 
tuned) systems it is the only available quantity. As the ef- 
fects of strong disorder are poorly understood at present, 
it is important to study very clean systems. The goal of 
this paper is therefore to establish a set of predictions for 
the resistivity within a SDW approach for weakly disor- 
dered metals close to an AFM QCP. The interplay of 
the relevant scattering mechanisms is studied within a 
semi-classical approach using a Boltzmann equation. The 
mod el and the Boltzm ann equation is defined in section 
II A . In section II B the analytic solution of the Boltz- 
mann equation is presented in the limit of low tempera- 
ture and weak, but finite disorder. The scaling properties 
close to the quantum critical point are emphasized and 
numerical solutions of the transport equations are used 
to analyze crossover regimes and non-universal behavior. 
Orbital effects of a magnetic field are studied in section 
III. Large non- linear effects can be used as a tool to inves- 
tigate the strongly anisotropic scattering from the spin 
waves. We will argue in section IV that these calculations 
easily explain some of the key observations in high purity 
single crystals close to an AFM QCP: the temperature 
dependence of the resistivity changes fijom,Ap ex T 15 to 
Ap oc T for cleaner and cleaner systemsDtm. We propose 
a number of stringent tests for our picture and discuss 
how the orbital effects of a magnetic field can be sepa- 
rated from spin-contributions. We conclude in section |v| 
by commenting on the relevance of this work to a wider 
class of problems. 



II. RESISTIVITY AT THE QCP 

A. Model 

Following HertztaEj a spin-density antiferromagnetic 
transition in a metal can be modeled by an effective 
Ginzburg-Landau theory defined by the action 
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where $ is the order parameter which fluctuates in space 
and (imaginary) time, /3 — 1/fcgT is the inverse tem- 
perature, bj n = 2im/ (3 are bosonic Matsubara frequen- 
cies. The fluctuations are strongest near the ordering 
vector Q of the AFM. The mass r controls the distance 
from the QCP. The Landau damping term linear in \uj\ 
is due to the effective scattering of quasi-particles from 
spin-fluctuations. A large number of particle-hole pairs 
is created by a spin-fluctuation if Q connects different 
parts of the Fermi surface (Q < 2fcp, see F-igJfl]). We will 
not discuss the singular case "Q = 2k p"i3r3 where the 
Fermi velocities Vk and Vk±Q are parallel or "Q > 2k p" 
where the spin-fluctuations do not couple directly to the 
quasi particles. 

Due to the ohmic Landau damping, the effective field 
theory (pi) is characterized by a dynamic exponent z = 2 
(w ~ q z \ an d is for d — 3 above its upper critical 
dimensionliSlia, d + z > 4. Therefore the dynamic suscep- 
tibility is mainly determined by the gaussian part of (jl]) 
and the spin-fluctuation jS^xj^trum in the paramagnetic 
phase can be modeled by 



XqM = X-qM 



1 



iM + ^-Wr' 



(2) 



where qo ~ kp and T are characteristic momentum and 
energy scales and £ is the correlation length of the spin 
fluctuations. The dispersion w q > vanishes at the or- 
dering wave vectors Qi and varies on the scale q$. For 
simplicity, we wilLmainly consider "isotropic" momen- 
tum dependences^ u> q 



((q ± Qi)/qo) 2 (appendix |D| 
defines the models used for numerical calculations more 
precisely). Moderate anisotropy influences our results 
only slightly, as will become clearer below. £ is the 
AFM correlation length which divergeaL§E3 at the QCP 
as l/£ 2 oc UT 3 / 2 . For the purposes of our numerical 
calculations we set 
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with c = 1. The T-dependence of £ does not affect the 
low-temperature properties (below Eqn. ( p8| ) a remark- 
able exception is discussed). We use c = 1 to model the 
(non-universal) destruction of the spin fluctuations at the 
temperature scale I\ The parameters in (0) can directly 
be obtained from inelastic neutron scattering. Typically, 
it is expected that in a heavy Fermion system, T is of the 
order of the coherence-temperature or Kondo energy Tjc- 
Close to the QCP, r is a linear function of the tuning pa- 
rameter, e.g. r oc p — p c in a pressure tuned experiment. 
Transport can be treated within a simple (quantum-) 
Boltzmann approach, because the spin-spin interactions 
are irrelevant in the renormalization group sense and be- 
cause furthermore the concept of Fermi quasi particles is 
still validEJ in 3 dimensions. For small and static elec- 
tric fields the transport equations take the usual form 



of fermions scattering from bosonic excitations. These 
equations are further simplified by the assumption that 
the spin-fluctuations stay in equilibrium, i.e. we neglect 
drag effects. This approximation implicitly assumes the 
presence sufficient momentum relaxation e.g. by strong 
Umklapp scattering. In the linear response regime, the 
quasi particle distribution /k = /£ — $k(9/£/<9ek) is lin- 
earized aroimd the Fermi distribution /£ an d the collision 
term readscll 
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Here gf and g| are transition rates for impurity scat- 
tering and inelastic scattering from spin fluctuations, re- 
spectively, and n° is the Bose function. 





p, x, B .... 



quantum critical point (QCP) 
FIG. 1. Left figure: In a metallic AFM the magnetic order 
is suppressed as a function of some control parameter which 
can e.g. be pressure, doping or magnetic field. Deviations 
from the usual Fermi liquid behavior show up close to the 
QCP. Right figure: Near the transition to an antiferromagnet 
with ordering vector Q, the scattering on the Fermi-surface 
is enhanced along "hot lines" connected by Q. This strong 
scattering equilibrates the distribution function (shown for an 
electric field E parallel to Q) in a region of width Afc (see the 
main text and appendix W for details) . 

In the following we will show that the interplay of the 
two scattering mechanisms is highly non-trivial because 
of their completely different momentum dependence. On 
the one hand, impurity scattering is very isotropic and 
is equally efficient on the whole Fermi surface. There- 
fore, the distribution function in the presence of a small 
electric field E is smooth <£>k °c VkE if impurity scatter- 
ing dominates (solid line in the right part of Fig. |lj). On 
the other hand the AFM susceptibility is strongly peaked 
around the ordering vectors Qi. Accordingly, a quasipar- 
ticle with an energy ek ~ p Wlu scatter efficiently from 
the spin fluctuations only near "hot lines" on the Fermi 
surface where £k H ~ £k H ±Qi ~ P- This strong scattering 
tends to equilibrate the distribution function at points 
Ilh and k# ± Q; in a region of width Afc (dashed line 
in Fig. |l|), where Afc depends on the relative strength of 



impurity- and spin-fluctuation scattering. The tempera- 
ture dependence of Afc and the distortion of the distri- 
bution function close to the hot lines in a magnetic field 
is the main origin of the anomalous transport properties 
which are discussed in this paper (see appendix [A^ for a 
simple qualitative calculation of Afc). 

For convenience, we will use in the following the dimen- 
sionless quantities t, x, r and b to measure the effective 
temperature, the amount of disorder, the distance from 
the quantum-critical point and the strength of the mag- 
netic field. 
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where pa = xpu is the residual resistivity and pu is 
a typical high-temperature (t « 1) resistivity which is 
defined below. One can approximately identify x with 
the inverse of the residual resistivity ratio (RRR). B 
is the typical magnetic field, which is necessary to see 
Shubnikov-de Haas oscillations at t = 1 
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The Fermi velocity vp is defined by an average over the 
Fermi surface u|/3 = 2/N F J dk(v k n) 2 / (vk(2ir) 3 ) where 
n is a unit vector in the direction of the electric field. 
J dk/uk = /<5(£k ~ /i)c? d k is an integral over the Fermi 
surface in d dimensions and Np — 2 J dk/ (vk(2ir) 3 ) the 
density of states. 



B. Boltzmann equation and scaling 

The linearized Boltzmann equations with the collision 
term (tm can be written in the following form 



ff dk' 
v k E + (v k x B)d k $k = // FW($k - *k') 77T-, 
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where an integration over directions perpendicular to 
the Fermi surface has already been performed and all 
k-vectors and integrations are restricted to the Fermi 
surface jf dk/v^ = J d d k6(e\ i — n) where d is the num- 
ber of dimensions. Currents are calculated from ji = 
J^$ k dk/(t; k (27r) 3 ). 

-Fkk' includes contribution both from elastic scattering 
from impurities and inelastic scattering from spin fluctu- 
ations with 
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I(y) is defined by equation 
totically given by I(y — > 0) : 
7r 2 /(3y 2 ). 
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For large temperatures and extremely large magnetic 
fields the solution of (0) depends on all details of the 
Fermi surface and only numerical solutions are possible. 
However, close to the quantum critical point the resis- 
tivity shows scaling behavior and depends only slightly 
on the geometry and the Fermi velocities along the hot 
lines. In the following we will derive an approximate an- 
alytic solution of the Boltzmann equation in this regime, 
the approximations used are exact in a well-defined limit 
of weak disorder, low temperature, weak coupling and 
small magnetic fields which will be specified below. Some 
of our qualitative results, especially in the absence of a 
magnetic field, can be derived by using a simple varia- 
tional approach (appendix H) or even in a simple relax- 
ation time approximation which we present in appendix 
0. The relaxation time approximation or, equivalently, 
the omission of vertex correction is completely wrong in 
certain geometries and leads to large errors especially for 
the Hall coefficient as explained in appendix |B| 

At T = the scattering is purely elastic and the 
quasiparticle distribution function is given by $ k (B) = 
$ k (S,T = 0) with $°(0) = eEv k /( 5i 2 mp A^). Our strat- 
egy is to keep track of the change of $k close to the hot 
lines, where the strong spin-fluctuations tend to equili- 
brate the distribution function at the points k and k±Q^ 
(Fig. [j]). For sufficiently small temperatures t 2 < x the 
distribution function will remain undisturbed in the "cold 
regions" far away from the hot lines where the inelastic 
scattering is of order t 2 and the elastic scattering domi- 
nates. It is therefore convenient to write the distribution 
function in the following forrncJ 
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For strong spin-scattering p k approaches a value close 
to 1/2 near the hot lines but vanishes further away. It is 
helpful to use a coordinate system, where vectors on the 
Fermi surface k = k# + kj^ are split up in a vector k# 
on the hot line and a perpendicular vector kj^ _L k# , v k . 
After a rescaling of the momenta using k = Kqo\/t, in the 
limit t < y/x the Boltzmann equation takes the following 
form 

t {d-i)/2 J dK ^ pK + PkH±Q+K , _ i)i(r/ t + AW) 

+yfuHC]t H xp K = c kff ^(v k x B)8kPk (10) 
with I(y) = J °° dzz 2 n a {z)(l + n a {z))/(y 2 + z 2 ), k' = 
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(2Tr) d N F ^/v kH v ]iH±Q /(TTq$), c / k[i = v kH v kH±Q /v 2 F andEl 

ag = w k +Q/u k . B and v k are unit- vectors in the direc- 
tion of B and v k , respectively. The strength of the mag- 
netic field is measured in units of a typical field Bq defined 
in (H) and (0) . One can think of b as the typical angle by 
which an electron is deflected within a typical scattering 
time tm = (2n) d VF / (2g 2 q ~ l ) . We are only considering 
magnetic fields with b <C x < 1 (i.e. uj c t <C 1), where 



Shubnikov-de-Haas oscillations are absent at low temper- 
atures. In (JHj) we have neglected sub-leading contribu- 
tions e.g. of the form p^d^^l, which are small in this 
regime because p k varies stronger with momentum than 
$° and all inelastic corrections come from small areas 
on the Fermi surface. This is not true in the ultra-clean 
limit x <C t 2 which has to be covered separately (e.g. in 
section [II A ). For completeness we discuss some of our 
result in d dimensions, however, we are mainly interested 
in three dimensions. Note that the Boltzmann equation 
approach presumably breaks down in d = 2 as the "hot" 
electrons acquire a lifctimoSil <x ,f/T, but it was neverthe- 
less widely used in literatureEIrEj. Qualitative features 
like the interplay of hot and cold regions might still be 
correctly described by our approach even in d = 2. 

From din) it is clear that only the combinations 
i( rf-1 )/ 2 /x, t/r and b/(xy/i) determine pk- Accordingly, 
the resistivity obeys a scaling relation 



(11) 



for t <§C \px < 1, r < 1 and sufficiently small b, where 
/ depends smoothly on the details of the Fermi surface 
near the hot spots and the directions of magnetic and 
electric fields, pm is defined in Eqn. (0). The relation 
(O) becomes exact in the scaling limit 
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The scaling limit serves as a convenient regime for an- 
alytic calculations. An experimental verification of the 
predicted scaling would be the most precise test of the 
underlying models. We will calculate the scaling func- 
tion / in the following paragraphs in detail. However, for 
som e app lications, e.g. the Hall effect discussed in sec- 



tion 



III A . the deviations from scaling are important and 



are calculated from a full numerical solution of the Boltz- 
mann equation ([7|) . A reader interested on ly in the qual- 
itative results can jump directly to section [V A, where a 
short overview of the asymptotic behavior of fa is given 
(the off-diagonal comp onents fi^j vanish in the scaling 



limit, see section III A where the Hall effect is discussed). 



The two dimensional coupled integro-differential equa- 
tions ( |l0| ) can be considerably simplified by realizing, 
that p k varies smoothly parallel to the hot lines but the 
kernel J(k[i) restricts the integrations to small values of 
k||. This allows us to replace p^ h ±q +K ' = Pk^+k' k' x by 
Pk«,k' m ( JIOJ ) and to perform the k',, integration. For the 
same reason we can neglect the contribution proportional 
to d K .,pK. All these approximations are valid in the scal- 
ing limit (O). It is therefore sufficient to solve a family 
of two coupled one-dimensional integral equations which 
depend parametrically on k# : 
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where 9b is the angle between B and a vector parallel to 
the hot lines at k# and 
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with 



r/t 



— 2kk' cos 0, where <f> ^ is 



the angle between v k/f and v kff -tQ. The equation for 
p K is obtained by replacing an by 1/ajj and exchanging 
p K and p K in (J13|). The boundary condition for (fill) are 
Pk,Pk — ► for k — > ±oo by construction. 
Currents are calculated from 
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where j = 2 J v k $ k d 3 k/(w k (27r) 3 ) is the current in the 
absence of inelastic scattering, §. dkn denotes the line in- 
tegral along the ith hot line and the expression is summed 
over all hot lines i. Note that p K is a function of k# . In 
the scaling limit t« i, the magnetic field dependence 
of $ k ex v kH E can be neglected in the second term in 
( |l6| ) . Using the symmetries of the Boltzmann equations 
one can derive the convenient expression 
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with Av l = vl H - vl H±Q . 



C. Transport for B - 

In the absence of a magnetic field, two different regimes 
emerge depending on whether ±he elastic scattering dom- 
inates close to the hot linecj'El If impurity scattering 
dominates (t < ruin [a;, \/rx\), it smears out the quasi 
particle distribution (Afc rs in Fig. nj). p K in (|13| ) is 
small and in leading order in p K one obtains 
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In this dirty limit, we recover the well-known T 5 resis- 
tivity at the QCP. This result is usually— derived by an 
average over the inelastic scattering ratetll. We want to 
stress that this approximation can be applied only in the 
presence of strong impurity scattering and gives wrong 
results for cleaner systems. 



nonuniversal regime 
t > sfx, p ~ t 2 




FIG. 2. In the scaling limit t,x,r — > 0, t/x,r/x —> const. 
the resistivity is "universal" Ap/pia = t 3 ' 2 f(t/x, r/x) for 
t < \fx,r < 1, where t oc T, x oc po and r oc p — p c 
measure the temperature, the amount of disorder and the 
distance from the QCP. The plot shows the qualitative be- 
havior of the scaling function f[a,/3] in the various regimes. 
For mm[x,^/rx < t < yfx] the resistivity rises linearly with 
temperature. Thermodynamic quantities show a crossover to 
Fermi liquid behavior at the scale t — r (straight line), while 
in transport a T behavior is recovered at a much smaller 
scale min[:r, y/rx\. The dashed line serves as a reminder that 
at lowest temperatures effects which are not included in our 
approach become impiestant, e.g. interference effects of dis- 
order and interactional ipia disorder induced change of the 
spin fluctuation spectrumEa. 

To derive (p"q) we considered for simplicity isotropic 
spin fluctuations. In the case of a moderate anisotropy 
our results have to be changed slightly. For example, 
if the anisotropy of the spin fluctuations is described 
by QigijQj/qQ, then the denominator of the integral in 
(|l|) has to be replaced by u kff « kH ±Q dct 1/2 JV kf g i3 \s?. 
where the k 1 is a unit vector parallel to the hot line and 
k 2 , k 3 are unit vectors perpendicular to k 1 and either 
v kif or v k/f ±Q. The determinant is always finite because 
we do not consider the case v k ||v k +q ("Q = 2kp"). 
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FIG. 3. Effective exponent of the resistivity at the QCP 
(r — 0), defined as the logarithmic derivative of Ap(T). At 
very low temperatures, the "dirty-limit" exponent 3/2 is re- 
covered. However, in the experimentally accessible low tem- 
perature regime smaller exponents are to be expected for 
rather clean system (x < 0.1). 

In the limit y/x > t > min[a;, y/rx], the resistivity rises 
linearly with temperature. The origin of this effect is 
that the spin fluctuations equilibrate a region of width 
Afc w t/y/x (Fig. [I]) as we will show in Eqn. ( |29|) (or 
appendix |A|). Therefore we obtain p ~ x(l + t/y/x) or 
more precisely 
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where the only slightly varying function S^ — 
J_ d,Kp K ,(a,a) « 7T is calculated from the solution of 
the two coupled one-dimensional integral equations (see 
Eqn. @) 
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{Pq + Pq' - 1)(1 ~ a 2 ) 

(q 2 /a + q' a — 2qq'a) 3 / 2 

{Pq + Pq' - 1)(1 ~ a') 

(q 2 a + q l2 /a-2qq , af/ 2 



= 







which in the limit a — > ±1 are solved hy p q = 1/((1 + 
a) + q 2 /2^c7) (see Fig. |). 

Omitting all prefactors we obtain in the scaling regime 
r< l,t< y/x 
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' d— 1 



(21) 



f t d / 2 ,r<t<a; 2 / d - 1 

^axji", y/xr^^] < t < y/x 
t 2 /r 2 ~"i , t < min[r, y/xr^^] 



2 4-d 



where for completeness we have given the result in arbi- 
trary dimensions d (see appendix [A] for a simple deriva- 
tion of these results). A summary of the scaling behavior 



in d = 3 is shown in Fig. |[ It is worthwhile to point 
out, that even away from the QCP, for r > 0, there is a 
large regime with a NFL resistivity (Fig. [| and ||). To 
actually observe the T 2 -term one has to consider very 
low temperatures t < \frx. This has to be compared to 
thermodynamic quantities, where the crossover to Fermi- 
liquid behavior is expected at a higher scale t ~ r. In 
a very clean system these scales can be quite different 
and anomalous transport should be observed in regimes 
where thermodynamic quantities display typical Fermi 
liquid properties. 
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FIG. 4. Effective exponent, defined as the logarithmic 
derivative of Ap(T), for a fixed distance from the QCP 
r — 0.01 and various values of disorder. In high purity sam- 
ples the crossover to an exponent 2 can be seen only at lowest 
temperature t <S min[:r, y/rx, r] jig). The bump at higher 
temperatures is a precursor of the effect that in ultra clean 



samples, x <Ct <C 1, we expect Ap ~ t 



In the scaling 



limit (113) the resistivity is linear in the intermediate regime 
y/rx < t < y/x (|2d). For the chosen parameters this exponent 
1 cannot be identified, however, a pronounced regime with an 
effective exponent less than 1.5 is sepjjh The inset illustrates 
the problem to extract the exponentou from the logarithmic 
derivative of Ap = p — po if po is not known. For = 0.02 
and x — 0.01 the logarithmic derivative of p — (1 — e)po for 
e = 2 — 10% (dashed lines) is compared to the "true" exponent 
(20) for e = (solid line). 



For higher temperatures t > y/x, larger resistivities 
Ap 3> po or further away from the QCP r ~ 1 the re- 
sistivity is much less universal and one has to rely on 
numerical solutions of the Boltzmann equation. In the 
ultra-clean limit x < t 2 -C 1, the distribution functioii 
<I>k is suppressed in a large region around the hot lincstll 
(Afc « O(kp) in Fig. m. Therefore only non-singular 
scattering dominates the resistivity and we obtain 
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(22) 



with a prefactor which to leading order is independent of 
the distance from the QCP r but depends on all details 
of the Fermi surface and of the interactions far from the 
hot lines. Hlubina and RiceEZI were the first to discuss 



this regime where quantum-critical effects cannot be seen 
in the transport because the long-living quasi particles 
in the "cold regions" short-circuit all contributions from 
the hot lines. Our numerical solutions, however, suggest 
that this regime is experimentally not accessible. Fig. H 
shows the logarithmic derivative of Ap, which defines an 
effective exponent a, Ap ex T a . Even for an extremely 
clean system with a residual resistivity ratio (RRR ~ 
1/x) of the order of 10 4 , there is never an extended range 
of temperatures, where this effective exponent is close to 
2. On the other hand we see that in moderately clean 
systems with x < 0.1 there is always a regime, where the 
effective exponent is smaller than 1.5. Sufficiently clean 
systems clearly show the linear resistivity predicted by 
( |20| ) over a large range. 

The full scaling function for r — and the deviations 
from scaling are shown in Fig. |ll|. 
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FIG. 5. Effective exponent, denned as the logarithmic 
derivative of Ap(T), for a relatively clean sample with 
x — 0.01 and various values of r oc p — p c . The resistiv- 
ity is calculated from a numerical solution of the Boltzmann 
equation for model A defined in appendix O. The inset shows 
p(t) for r = 0,0.1,1. 



Fermi surface where the scattering is relatively strong 
(or weak) then non-linear effects will become important 
for uj c t > Ak/kp. Close the AFM quantum critical point 
the width of hot lines Afc vanishes in the low-temperature 
limit and the quasi particles can be driven by a moderate 
magnetic field over the hot lines. This will lead to large 
effects in the magnetotransport. The magnetotransport 
in nearly antiferromagnetic two dimensional metals has 
been studied beforeSEIZl in detail in the context of high- 
temperature superconductors. However, the authors did 
not focus on the role of weak disorder, which we find 
at least in three dimensions to be essential to describe 
the low-temperature the magnetotransport both quanti- 
tatively and qualitatively. 
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FIG. 6. Hall effect for model A, calculated from a numerical 
solution of m Extremely small amounts of disorder strongly 
influence the Hall effect. Very clean samples are necessary to 
observe this effect. Both the sign and the size of the effect 
depend on details of the band structure. 



A. Hall effect 



III. MAGNETOTRANSPORT 

The characteristic feature of the spin-density wave sce- 
nario is, that only the quasi particles close to "hot lines" 
are strongly affected by the spin-fluctuations. This has 
to be contrasted with a situation where the Kondo effect 
is destroyed close to the QCP which would affect the full 
Fermi surface. It is therefore important to have a probe 
for testing the presence of sharp structures on the Fermi 
surface. 

The (nonlinear) magnetoresistivity can be used to in- 
vestigate precisely this question. By a magnetic field, 
the quasi particles are driven parallel to the Fermi sur- 
face by the force vtxB. Generally, non-linear effects in 
the magnetotransport show up if the quasi particles are 
able to circle the Fermi surface without being scattered, 
i.e. for lo c t > 1, where uj c is the cyclotron frequency. 
However, if there is a small region of width Afc on the 



The Hall effect is extremely sensitive to small amounts 
of disorder (x < t 2 ) but essentially constant for t <C %/x. 

This can be seen most easily from the symmetry prop- 
erties of the Boltzmann equation (13) in the scaling limit. 
The left-hand side of ( |13| ) is even in n, while the right 
hand side is odd and linear in the magnetic field. There- 
fore J_ p K is even in B and there is no contribution to 
the Hall effect from the hot lines in the scaling limit ( p^ ) 
where ( |l3| ) and (|]J]) are valid. This does not imply that 
the Hall effect is unaffected by the singular scattering 
close to the QCP. It is on the contrary an extremely sen- 
sitive probe to very small amounts of impurity scattering. 
The Hall constant at T — is different in the ultra-clean 
and the dirty limit Rh {t — > 0, x < t 2 ) ^ R H {t — > 0, x = 
const.). In the former limit the elastic impurity scat- 
tering dominates in the latter the magnetic scattering. 
The quasiparticle distribution $k has a completely dif- 
ferent momentum dependence in these two limits. One 



consequence is that a different average over the effective 
masses enter the Hall constant. Small amounts of dis- 
order x ~ t 2 change the distribution <I> k completely and 
give rise to the strong features in the Hall effect shown 
in Fig. pi The calculation is based on a numerical solu- 
tion of the full Boltzmann equations (|7|) for small mag- 
netic fields. Even disorder at the level of e.g. x = 1CP 4 
leads to a strong structure in the Hall effect at quite high 
temperatures. Extremely clean samples are necessary to 
observe the strong temperature dependence of Rr- Both 
the sign and the size of these effects are non-universal 
and depend on the details of the band structure. 



B. Magnetoresistivity — small fields 

While all corrections to the Hall effect cancel in the 
scaling limit (13), the magnetoresistivity is strongly in- 
fluenced by the hot lines and can be used to investigate 
the quantum-critical properties in more detail. 

We first discuss the dirty limit x > t (d — 3h p K 
is small in this case allowing an expansion of ( |13| ) in 
p K . In order 6 2 , the distribution function changes by 
Ap K sa (6cose B c / kff /(x v / tc kif% /ai")) 2 3^ where p?is 
the distribution function for 6 = 0. The resistivity ( |17| ) 
is proportional to J p K , but J Ap K vanishes obviously. 
Therefore it is useful to integrate (Oh 



dnp K 



dndK'G KK > (1-Pk-Pk 1 )- 



(23) 



We calculate the leading correction to / p K from the right 
side of (B3), using that p° w £/(x^/a_ffC kif ) J G KK 'dn': 



A p K 



1 + a H ( bcosQ BC / kH 



d K (d K p° K ) 2 (24) 



and obtain for r < t < min[a;, \/rx\, 6 <C x^Jr 
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(25) 



where B is a unit vector in the direction of the magnetic 
field. It is surprising that the leading low temperature 
correction in this Fermi liquid regime t <C r is not pro- 
portional to T 2 but starts with T A . This is due to the 
above mentioned cancellations and is actually valid only 
for t > T\/x~ as shown below. 

There are a number of other contributions to the mag- 
netoresistivity which vanish in the scaling limit but might 
be dominating in the experimentally relevant regime. 
One important correction is due to the suppression of the 
AFM in a magnetic fields, we will discuss these effects in 



detail in section IV C. Before analyzing the magnetore- 



sistivity it is usefull to subtract the B-dependence of the 
resual resistivity Ap ~ b 2 /x which is the largest correc- 
tion to our results. Therefore we use in this section the 
definition 



Ap = 09(6, t) - p(0, t)) - (p(b, 0) - P (0, 0)) 



(26) 



Other sub-leading contributions become important for 
higher magnetic fields or in the case of Eqn. fl25| ) 
for low temperatures. They are missing in ( |10| ) be- 
cause we neglected terms of the form dk(^ H — $k H ±Q)- 
We estimate these corrections in the regime uj c t <C 1 
and Ap < p using Kohler's rule (p(T, B) — p(0, B) oc 
B 2 Ap(T,0)/p(0,0) 2 ). In the regime discussed in ©, 
t hey are of the order of Ap oc b 2 t 2 / (tJtx 2 ) and therefore 
(25) is valid only for t > t* = r^/x. In the scaling limit 
(12), t* vanishes. 

The behavior in the disorder dominated quantum- 
critical regime t < x, r -C t is more complicated. 
For small momenta, G KK i is approximated by G KK i ~ 

l/(2-i r/t + k 2 + k' 2 — 2kk' coscj)), where <f> is the angle 
between v kff and v k „ + Q. In leading order in t/x for 
6 = and k < 1 Eqn. J7] is solved by 

p° K ~t/x f G KK ,dK f ~ t/(2x)ln[l/(r/t + « 2 )]. (27) 



where the prefactors and <f> dependence in the argu- 
ment of the log have been omitted. The approxima- 
tions leading to (E7p are valid for p° K <C 1, therefore they 
break down in an exponentially small regime r/t + k 2 < 
ce~ x K 2t \ where c is a constant of order 1. 

With all prefactors, which are calculated using (|24|) , 
for e - ^ 2 *' <C r < t < x, the inelastic contribution to 
the magnetoresistivity is given by 
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(28) 



where the geometrical factors are accidentally exactly the 



same as in 



. It is quite surprising, that even in the 
limit r<f the quadratic magnetoresistivity diverges pro- 
portional to 1/^/r. What happens if a system is directly 
at the QCP, i.e. in the case r — 0? In the scaling limit 
©, the cutoff r/t in @ is replaced by e~ x/{2t) and 
the magnetoresistivity is proportional to 6 2 i 1 ' 5 e a;/ '( 4 *) /x 3 . 
This exponentially large correction is probably not exper- 
imentally observable. Either higher-order effects in the 
magnetic field will become relevant (to be discussed in the 
next section) or one has to take into account the tempex^. 
ature dependence of r. The low-temperature correctionES 
Ar rs c'i 3 / 2 can formally be neglected in the scaling limit 
and usually leads only to small sub-leading corrections of 
the order of Ar ft m c'\/t both in transport and thermo- 
dynamic quantities. The prefactor c' is proportional to 



the spin-spin interaction U (|l|). The magnetoresistivity 
at the QCP in the dirty limit is a remarkable exception 
of this rule: for r <C t <C x and c' ' \fl > e~ x '( 2t ' the tem- 
perature dependence of r determines the B 2 correction 
to the resistivity and we expect Ap ~ b 2 t 125 /(Vdx 3 ) 
( Pq ) ( see Fig. [7j). An experimental confirmation of this 
prediction would be an interesting way to investigate the 
temperature dependence of l/£ 2 , a quantity which can 
otherwise be measured only by high-precision inelastic 
neutron scattering. 
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FIG. 7. Log-log plot of the temperature dependence of 
the magnetoresistivity Ap/b 2 — Apx/b 2 at the QCP in a 
small magnetic field b for model A, calculated from a nu- 
merical solution of (|7j). The "Kohler's rule" -contribution 
Ap K - Apit = 0, b)(p(t — 0,b- 0)/p(t, b - 0)) has been sub- 
tracted. In the dirty limit, t < x, one can indentify the correc- 
tion @ with Ap ~ b 2 t 2 /(x 3 ^r) ~ fo 2 t 5/4 /x 3 , r w i 3/2 at the 
QCP. The magnetoresistivity in the clean regime, x < t < yjx 
is much larger with Ap ~ b 2 /(y/xt) (see Eqn. (p3)). In the 
ultra-clean limit, x < £ 2 , we obtain Ap ~ 1/t 2 and to leading 
order no effects of the QCP can be seen. 



The solution of the Boltzmann equation ( |13[ ) in the 
regime max [a;, ^/rx] « ( < y/x is dominated by large 
momenta k,k' 3> max[l, \/r/t]. We therefore use the 
asymptotic form of G KK > (15h and rescale the argument of 

p K by k = qyTTt/(6xck H sin 2 <j>)/y/aH and the argument 



of p K > by k' = q'yirt/iSxCkj. 
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{q 2 /a H 

dq'(p q 



q' 2 a H 



2qq' cos(j)) 3 / 2 
Pqi - 1) sin 2 



tva; 



(g 2 aif + q' 2 /a H - 2g</ cos 0) 3 / 2 ^v^ 



= b -^d q p q (29) 



with £k ff = cjj. | sin^lv^rJcosOs/^/Trcka • For small mag- 
netic fields, p q decays on the scale 1 (both an and <p & r e 
finite for a generic Fermi surface). This justifies our ap- 
proximation to expand G KK i for large momenta in the 
regime \Jtjx > max[y / r/i, 1], i.e. for t > max[y/rx,x]. 
In the original units the width of pk is approximately 



Afc rs qot/y/x and by scaling we obtain the linear resis- 
tivity shown in Eqn.(|20j). 




FIG. 8. Shape of p q in the regime b = 0, 
max.[y/xr, x] < t < yfx for various values of olh — «k±q/vk 
and cos<?!> = v k ±q ■ v k /(vk±QWk)- The integral S — J p q ~ ty 
depends very weakly on the parameters. p q with k oc qt/y/x 
describes how strongspin fluctuations equilibrate the distri- 
bution function $k (PI) an d *&k±q in a region of width t/y/x 
around a hot line. 



While the kernel of ( |29| ) is highly singular, the so- 
lution of the integral equation for b = is neverthe- 
less very smooth: for large momenta p q decays propor- 
tional to \/q 2 and is finite and smooth for q — > with 
Po + Po = 1 to cancel the divergence of the kernel (see 
figure @). The perturbation theory in b/(y/xt)) is an- 
alytic and well defined (all eigenvalues of the matrix 
5(q — q')(l + f Gqq'dq') + G qq > are larger than 1) and 
we obtain for b < y/xt, y/x > t > maxf-y/rx, a;] : 
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(30) 



Although it is difficult to calculate the precise prefactor 
analytically, its structure is similar to the prefactor in 
(pq). In Fig. m the numerical solution of the Boltzmann 
equation shows the crossover from the dirty limit d2q) to 
the clean limit (|30"|). In the ultra-clean limit, x <C^the 
distribution function varies on the scale kp and Kohler's 
rule is valid in leading order and therefore we expect 
Ap~b 2 /t 2 . 



C. Magnetoresistivity — large fields 

The diverging prefactors of the B 2 -corrections already 
indicate that any finite magnetic field will strongly influ- 
ence the resistivity when the hot lines get sharper and 
sharper while B drives the quasi-particles over the Fermi 
surface. 

For large fields, the main effect of the bd^pk-tevm in 
the Boltzmann equation (113) is to keep d^pk small. Due 
to the boundary conditions Pk—f±oo ~* and due to the 



fact that the kernel Gkk> decays for large momenta, the 
amplitude of pk decreases for large magnetic fields and 
Pk is broadened. This corresponds to the fact that mag- 
netic field drives the quasi particles over the Fermi sur- 
face hence smoothing the quasiparticle distribution. In 
a strong enough magnetic field we can therefore use per- 
turbation theory in pk and only the first term on the 
right-hand side of Eqn. ( |23| ) survives. In the scaling 
regime, the resistivity saturates in a large (orbital) mag- 
netic field at a value which is precisely given by the dirty- 
limit formula (lq). In this respect, a magnetic field and 
disorder have very similar effects, both smooth out the 
quasi-particle distribution. 

As the amplitude of pk decreases for large b, the 
main effect of a large magnetic field in the dirty limit 
t < max[i, \/rx\ is to suppress sub-leading corrections to 
the conductivity. The maximal change of J p K is calcu- 
lated from (ES 



A p K (b -*oo) 



t 



o i ^o 



dndK'G KK '(p K +p 



= (l + a H ) f°°dK(p°f (31) 

J — oo 

where p° K — t/(x^/aHCk H ) J G KK idn' is the distribution 
function in the dirty limit in leading order in t/x. 

The effect of an intermediate field can be calculated by 
solving the Boltzmann equation in a finite field in pertur- 
bation theory in t/x by neglecting again the G KK > (p K + 
p K >) terms. 



Pk 



dn'p%'- 



= -(fc'-fc)/«i, 



Kb 



(32) 



with Kb = b/(xVi)(c[ i /(i/cwckjj)). The main effect of 
the magnetic field is to smear the distribution function 
on the scale Kb- 

In the Fermi liquid regime t < min[r, ^/rx\ the scale of 
the distribution function p a K ~ (t/x)/(r/t + k 2 ) is given 
by \/r/t, accordingly the small-field formula (Bq) can be 
applied for Kb < yr /t or b < \prx. As discussed above, 
in a large magnetic field Ap defined in (J26|) saturates at 



Ap(b > y/?X, t) 



t 4 



Pm 



r 3/2- 



(33) 



and the weak-field expansion (|2jj) connects smoothly to 
this regime. As above, we can estimate the range of valid- 
ity of ( p3| ) by comparing it to neglected corrections of the 
order of Ap - b 2 Ap(0,t)/p 2 ~ b 2 t 2 /{■s/rx 2 ). Therefore 
( p3| ) is valid for b < ty^x/r, t < min[a;, y/rx\. 

The dirty quantum-critical regime r < t < x has a 
richer structure, as there are two scales in the distribu- 
tion function p Q K (pTT): 1 and max[^/r/i, e~ x ^ 2t ^]. For 
1 > Kb > max[yV/i, e^ x ^ 4t ^] the ma gnet ic field, or more 
precisely k" 2 , takes the role of r/t in ( |28|) and provides a 



cutoff to the logarithmic divergence of p° (pPj 
fore expect for r < t < x, max[x-\/ie~ x '( 2 ^~. 

x\fi 



We there- 
frx\ < b < 



MM) ,t 2 
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(34) 



and the magnetoresistivity will saturate for b > xy/i at 



A P (b,t) t 5 / 2 
Pm x 



(35) 



This has to be compared to the largest sub-leading cor- 
rection of the form b 2 t 15 jx 2 (estimated from Kohler's 
rule). Therefore equation d33) is valid for b < \/tx. 

All the above discussed effects in the dirty limit are 
very small on an absolute scale because the disorder has 
already smeared out most of the features of the distri- 
bution function and can only be observed because it is 
easy to measure the difference p(T, B) — p{T) with high 
precision. Much stronger effects are to be expected in 
the cleaner regime max[:r, ^/rx\ < t < t/x, where the 
zero-field resistivity is linear in temperature. 



Q. 




FIG. 9. Schematic plot of the temperature dependence of 
the resistivity. In an ultra-clean system, the resistivity would 
rise proportional to t (thin dotted line). In a system with a 
small amount x of disorder the resistivity at the QCP (thick 
solid line) at zero magnetic field rises with t 3 ' 2 at very low 
temperatures t < x, but shows a linear behavior Ap ~ ^fxt 
for x < t < yfx. If the system is tuned to the quantum-critical 
point in a finite field, a behavior similar to the dot-dashed line 
is expected: Ap oc t 3 ' 2 should be observed over a larger range 
t < b 2 ' 3 . In an intermediate regime b 2 ' 3 < t < b/yfx the 
resistivity depends weakly on temperature (the thin dashed 
line corresponds to Ap oc t 3 ' 2 ). Note that a different behavior 
is expected if a finite magnetic field is applied to a compound 
which is quantum critical for b = as b drives the system 
away from the QCP. 

For b < ^Jxt, the effect of the magnetic field can be 
treated in perturbation theory and ( |30| ) is valid. To an- 
alyze the regime \fxt < b it is useful to rescale the mo- 
menta in (E3) using q = —bc] iH /(ty/x)ky?a>H so that the 
Boltzmann equation reads 



-d k Pk 




-2kk'cos(j)f/ 2 
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(36) 

The corresponding equation for §y is obtained by re- 
placing an by 1/a.ff. The kernel of the integral equa- 
tion is strongly divergent J Gkk'dk' oc 1/fc 2 . The con- 
tribution from the integral is, however, small for large 
magnetic fields and sufficiently large momenta k > fco 
(a precise estimate is given below) and (p6| ) is solved by 
Pk^o — c±e~ k ^° r " . For k < the prefactor c_ has to 
vanish due to the boundary condition p\k\->oo — 0. For 
k — > the divergence of the kernel dominates the inte- 
gral equation and forces the solution to take the form 
Pk~>o + Pk'^o — 1- For an = 1 one therefore expectsE3 
Pk^a = Pk'^a = 1/2- Generally, pk^o is a number close 
to 0.5 depending slightly on an,4> and b. For a large 
magnetic field, the rise of pk close to k — is indepen- 
dent of a.H andpfc^o(fr 3> \/xt) = Pk^o{b ^> y/xt) = 1/2. 
The prefactor of the exponential decay c+ for k > fco is 
therefore 1/2 and J Pkdk = 1/(2^/oh) or, after rescal- 
ing, Jp K dK = 6/(a;'v / t)|cos0 B |(c( CH /(2aHC kfr )). In this 
regime y^i < b < min[£ 3 ' 2 , t 2 /y/r], max[x, yfrx\ < t < 
y/x, the rise of the resistivity p— po is linear in |B| and in- 
dependent of temperature and disorder (the precise range 
of validity is calculated below) . 
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For higher magnetic fields the formula given above will 
break down, because (E9h is valid only for large momenta 
K > max[l, y/r/t] or, equivalently, (pq) is valid for k > 
m&x[x\/t/b, y/rxjb]. To estimate the range of validity of 
(pTJ), we solve ( pq ) more precisely for k < in a regime, 
where pk is still small: 
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Pk is large for k < fco ~ t 2 x/b 2 and Eqn. ( p?\ ) is 
therefore valid for fco > m&x[y/tx/b, y/rx/b] or b < 
min[£ 3 / 2 , t 2 /y/r\. For larger values of b the kernel is less 
singular and pk stays small for all momenta. Accordingly, 
the resistivity is dominated by the first term on the right- 
hand side of Eqn. (|23| ) and the magnetoresistivity satu- 
rates at a value, which is given by the dirty-limit formula 
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This expression is valid for min[£ 3 / 2 , t 2 /y/r\ < b < 
jmii[t 1 ' i x s ' i , \/ix 3 ' i /r 1 ' i ], where the upper limit esti- 
mates the regime where non-universal contributions fol- 
lowing Kohler's rule dominate. In Fig. |j the temperature 
dependence of the large- field magnetoresistivity is shown 
schematically. 



IV. DISCUSSION 

A. Overview over magnetoresistivity 

Before we discuss the experimental relevance of our 
calculations in detail, we collect the results for the mag- 
netoresistivity in the scaling limit (12). Most of the pref- 
actors, which are omitted here, can be found in the previ- 
ous sections. Only orbital effects of the magnetic field are 
taken into account, a discussion of the interplay of spin 
and orbital effects is given in section IV C . t and b defined 
in (H) are the dimensionless temperature and magnetic 
field, respectively, x is proportional to the strength of 



disorder and can crudely be identified with the inverse of 
the residual resistivity ratio and r measures the distance 
from the QCP. 

In the disorder dominated regime the absolute size of 
effects of an orbital magnetic field are quite small. Close 
to the QCP, for r <C t < x, the prefactor of the 6 2 -term 
is, however, extremely large and for very small magnetic 
fields a crossover to a linear field dependence of the mag- 
netoresistivity is predicted. For Ap = p(6, t) — p(6, 0) we 
obtain 



Ap 

PM 



+3/2 



t 5/2 



b 2 t 2 



X 

^5/2 

X 



,b< g x 



r.t 



,9x,r,t <b<xVt (40) 
,Xy/t < b < \/tx 



with g x , r .t ~ max[xy/r, x\fie~ x > ^'] . For very small g x , r ,t 
this expression can be used to measure the temperature 
dependence of the correlation length £ using r oc l/£ 2 . 

In the disorder-dominated Fermi liquid regime t < 
min[r, y/rx] < 1, the temperature dependent part of the 
magnetoresistivity is still enhanced, but it is probably 
very difficult to extract the t A contribution experimen- 
tally. 
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PM 



XV 
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b 2 t A 
i 4 



xr 



3/2 



, b < yfrX,t > Ty/x 



, yfrx < b < t 



(41) 



The largest effects can be observed for rather clean 
systems with max[x, y/rx\ < t < y/x, which show a linear 
resistivity for 6 = 0: 
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PM 



\fxt 



\fxt 



,b < y/xt 



, y/xt < b < min[i 3/2 , —=] 



(42) 



i 2 t 2 

iin[£ 3/2 , -=] ,min[£ 3 / 2 , —=}< b < 

x^minlt 1 ^,^] 



B. Experimental situation for B — 

According to our calculation, in d — 3, spin-fluctuation 
theory predicts that the resistivity in dirty systems rises 
with T 3 ' 2 , while in cleaner compounds a large regime 
with a linear resistivity is expected. This behavior is ac- 
tually seen in a large number, but not all, AFM quantum- 
critical systems. A resistivity proportional to T 1 ' 5 has. 
been repotted in CeCu^y^rCeNiGaJ, C<»Cu 6 - x AgjB 
CeNi 2 Ge 2 la, CePd 2 Si 2 , CelngB or CeCu 5 AiM All these 
systems are dirty in the following sense: in the regime, 
where the exponent 1.5 has been observed, the rise of the 
resistivity Ap is small compared to the residual resistiv- 
ity. 
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FI£J. 10. Resistivity as a function of T 1 ' 2 in CePd2Si2 taken 
fromQ (left figure) compared to our calculation within model 
A defined in appendix Ofor x = 0.01 (right figure). The insets 
show the corresponding logarithmic derivative of p(T) — p(0). 
If one compares the insets one has to take into account that 
the low-temperature behavior of the logarithmic derivative 
of Ap strongly depends on the exact determination of po (see 
inset of Fig. ]a). Below ~ 400mK CePd2Si2 is superconducting 
(lower inset). Note the offset of the line T = in both plots. 

For a small number of the above mentioned com- 
pounds, cleaner samples exist. For these high-purity sin- 
gle crystals, the resistivity seems to rise with an exponent 
smaller than 1.5. The cleaner the samples are, the closer 
the exponent is to 1. This qualitative trend coincides 
with the predictions of our calculations, jjEig. 0). Ex- 
amples for this are CePd 2 Si 2 , CeNi 2 Ge 2 ErEHj where sam- 
ples with RRR of the order of 100 exist. In clean sam- 
ples of Celn3, exponents smaller than 1.5 are observed 
over some temperature range, but at lowest temperatures 



the exponent crosses over to 1.5, as expected from the- 
ory. In Fig. [L0| the resistivity of a quite clean sample of 
CePd 2 Si 2 H is compared to a typical solution of the Boltz- 
mann equation. It is important to point out, that not 
only a similar effective exponent shows up in both the- 
ory and experiment, but that it is observed over a similar 
range O.lpo < Ap(T) < 10po- Only the effective scale 
r/f/Q seems to be approximately a factor 10 too large 
compared to Tx/kp, but a quantitative comparison of 
these high-energy scales is certainly very dubious within 
our perturbative treatment of a strongly correlated sys- 
tem. Furthermore, a quantitative comparison would re- 
quire a more realistic model based on some knowledge 
of the band structure and the actual width of the spin- 
fluctuation. 

The theory, presented in this paper, is certainly not 
applicable to strongly disordered systems, where the na- 
ture of the quantum-critical transition is changed due 
to disorder. It does not explain the linear tempera- 
ture dependence of the resistivity in CeCug-icAua; due 
to two reasons. The first is, that CeCus-zAu^ is rela- 
tively dirty at the QCP (x — 0.1). More important is 
the unusual nature pOf the spin-fluctuations, which are 
strongly anisptropicL3, and show an unusual frequency 
dependence^. 
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FIG. 11. Scaling plot of f(t/x,Q) = Ap(t/x)/t 3/2 at the 
QCP (r = 0) using a numerical solution of the Boltzmann 
equations for various amounts of disorder x. The straight 
solid lines show the limiting cases Ap oc t 3 ' 2 for i < 1 and 
Ap oc t for 1 < t <C y/x. Deviations from scaling are large 
for t/x > 1/y/x. 

The most precise test of the spin-fluctuation model 
would be to verify (or falsify) our scaling prediction (|ll| ) 
by comparing several very clean quantum-critical sam- 
ples with varying residual resistivity ratios. An exam- 
ple for the expected deviations from scaling is shown in 
Fig. H 

A somewhat unexpected result of our calculation is 
that in clean systems away from the QCP, a Fermi liq- 
uid form of the resistivity can be observed only at very 
low temperatures t < \frx. This has to be compared to 
thermodynamic quantities, where FL behavior sets in at 
a considerably higher scale t < r. In a clean system with 



12 



a RRR over 100, the crossover scale in transport quan- 
tities can easily be e.g. a factor 5 lower than in ther- 
modynamic quantities. This might explain claims based 
on transport measurements that the non- Fermi liquid be- 
havior in these clean-systems extends over a finite region 
away from the QCFm 



C. Spin- and orbital effects of B 

A magnetic field influences the resistivity near a quan- 
tum critical point in two quite independent ways: by 
spin- and orbital effects. The spin contribution typically 
suppresses the antiferromagnetic order. In the context of 
this paper, we will call these phenomena "spin effects" 
despite the fact, that in most cases orbital degrees of 
freedom are involved. On the non-magnetic side of the 
phase diagram and for B fields, which are not too large, 
the leading order spin-effects can supposedly be described 
by the magnetic field dependence of r, the variable which 
describes the distance to the QCP. For small fields, one 
expects in an AFM metal above the upper critical dimen- 
sion Ar ex B 2 . A controlled study of the effects of a finite 
magnetic field close to the AFM QCP in a metal does not 
exist to our knowledge, however, Ioffe and Milliscla have 
carefully investigated the uniform susceptibility without 
finding unexpected anomalies. A large magnetic field can 
induce other (first-order) transitions or change the nature 
of the magnetic fluctuations completely. 

The focus of this paper are "orbital effects" . We con- 
sider a regime where interference effects don't play any 
role and have restricted our attention to the classical 
v x B force on the electrons. This approximation is valid 
in the scaling limit (|1J), where uj c t -C 1 and the disor- 
der is so weak that weak-lecalization corrections set in at 
much lower temperaturesE3. 

In the scaling limit (Il2), i.e. for clean enough systems 
at very low temperatures in small magnetic fields close 
to the QCP, spin effects can be neglected. For b ex x 3 ! 2 
the correction Ar oc b 2 is small compared to x for suf- 
ficiently small x. In realistic heavy Fermion systems, it 
is, however, probably not possible to produce a sample 
which is clean enough so that only orbital effects show 
up in the interesting regime. The reason is that on the 
one hand the Lorentz force v x B is suppressed due to 
the large mass m* of the heavy quasiparticles and on 
the other hand spin-effects are enhanced due to the large 
susceptibility x ~ 1/^if oc m* /mo- 

A reliable estimate of these effects is difficult as the de- 
tails of the interplay of the Kondo [effect and the RKKY 
interaction are poorly understoodc2l For a crude esti- 
mate, we assume that all high-energy scales, including 
scattering rates and renormalized coupling constants can 
be approximated by Tr . Orbital effects are important for 
b > x 3 ! 2 or hbB > x 3 / 2 kBTKm* /mo and spin effects can 
be neglected for Ar < x or iibB < ^JxIubTk y/mo/m* . 
If this is qualitatively correct, only an extremely clean 



sample with x < (mo/m*) 3 / 2 , i.e. with a RRR larger 
than (m*/mo) 3 / 2 , shows orbital effects without the in- 
fluence of spin effects. This implies that spin-effects are 
present in most cases, because m* /mo is quite large in 
those clean stoichiometric systems where the quantum 
critical point is reached by pressure. 

Nevertheless, the orbital effects have a very clear signa- 
ture even in the presence of large spin-effects, especially 
in suitably designed experiments. The necessary mag- 
netic fields are not very large. In our units, Shubnikov-de 
Haas oscillations are expected for b > x at low temper- 
atures while strong effects in the magnetoresistivity in 
clean samples set in for b > x 3 ^ 2 , i.e. for magnetic fields 
which can easily be an order of magnitude smaller. 

The Hall constant is obviously not influenced by 0(B 2 ) 
corrections due to the suppression of the AFM and so Rh 
is not affected by the above mentioned problems. There- 
fore it should be very interesting to investigate the strong 
samp le dependence of the Hall effect predicted in section 
[II A. Unfortunately, only extremely clean samples-will 
show a strong temperature dependence (see Fig. ^)cil. 

For a detailed study of the magnetoresistivity in clean 
systems, it is necessary to minimize spin effects. In 
a system like CePd2Si2 the effective distance from the 
QCP can be tuned by pressure. Therefore, one can se- 
lect various pressures on the antiferromagnetic side of 
the phase diagram close to the QCP and suppress the 
remaining magnetic order by a suitable magnetic field. 
For such a set-up the spin effects are compensated by 
a change of pressure to a large extent especially in sys- 
tems with Ising symmetry. For not too strong fields the 
field-tuned quantum-critical transition should stay sec- 
ond order and one can compare several quantum-critical 
systems (r — 0) with the same amount of disorder x but 
different magnetic fields. Our calculation predicts, that 
in a magnetic field b > x 3 ^ 2 , the low- temperature resis- 
tivity rises proportional to t 15 . A schematic plot of the 
expected behavior is shown in Fig. |[ Ideally, one can test 
the scaling predictions by comparing samples of different 
quality. In an anisotropic system one should choose a 
direction of the magnetic field perpendicular to the easy 
axis and to one of the ordering vectors Q, because only 
the component of B parallel to the hot lines enters the 
scaling form of the resistivity. (This effect is described 
by the (vk,, x Vk,,±Q )B factors in our results, e.g. in 

A simple way to distinguish spin- and orbit effect is 
the sign of the magnetoresistivity. A suppression of the 
AFM fluctuations on the paramagnetic side of the phase 
diagram will reduce the amount of scattering and the 
resistivity will go down. Due to the orbital effects dis- 
cussed in the previous section resistivity increases be- 
cause b smears out the quasi particle distribution, there- 
fore minimizing the effect that cold regions short-circuit 
the hot spots. This can lead to a sign change of Ap = 
(p(B, T) - p(0, T)) - (p(B, 0) - p(0, 0)) either as a func- 
tion of disorder or strength of the magnetic field. Such an 
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observation could help to separate the two effects quan- 
titatively. 

The interplay of the spin- and orbit effects can also 
lead to less trivial effects. In the dirty limit close to 
the QCP e.g. the magnetoresistivity for small magnetic 
fields (H) is proportional to 1/s/r, Ap ~ b 2 t 2 /(x 3 ^/r). If 
spin effects suppress the AFM fluctuations with r « cb 2 , 
one can expect a regime where the magnetoresistivity 
rises linearly with B for rather small fields (cb 2 < r(b = 
0, tj), Ap oc bt 2 /(^/cx 3 ). It is probably not very easy to 
distinguish this effect from the linear b dependence due 
to orbital effects discussed in ( |34| ) with Ap ~ bt 2 /x 2 for 
b>xy/r(b = 0). 

ti high purity single crystals of CeNi2Ge2 Grosche et 
have reported measurements of the resistivity in mag- 
netic fields of a few Tesla. As explained above, a di- 
rect comparison of their findings with our calculations 
is difficult due to the interplay of orbital and spin ef- 
fects. However, an important qualitative trend of their 
results can easily be explained by our approach. In fi- 
nite magnetic field, they observe a clear T 2 dependence 
below a temperature T*(B). Obviously, the magnetic 
field drives the system away from the quantum-critical 
point. Remarkably, it seems that the magnetic field is 
much more effective in this respect than pressure, which 
also suppresses the magnetic order. For samples tuned 
by pressure towards the paramagnetic side of the phase 
diagram, a clear T 2 -rise of the resistivity could not be ob- 
served. This qualitative effect can be understood within 
our approach. In the case of pressure tuning, we expect 
Fermi liquid in a very clean sample only for very low 
temperatures t < \prx. In a sufficiently strong magnetic 
field, the crossover to T 2 should be visible at a much 
higher temperature t < r l ^\/b in agreement with the 
experimental trends. If we assume that r ex b 2 and use 
the above discussed crude estimates of the prefactors, 
we expect the crossover temperature T* to rise linear in 
B with ksT* ps psB^o/m*) 1 ' 4 . I want to emphasize 
that this estimate is unreliable and should be replaced by 
a suitable strong-coupling theory of heavy-Eermion sys- 
tems. Nevertheless, the experimental resultsQ seem to be 
consistent with this simple-minded estimate of scales. 



V. CONCLUSION 

In all situations where two scattering mechanisms com- 
pete, one anisotropic and affecting only parts of the Fermi 
surface, the other isotropic, the transport is quite subtle 
and the interplay of the two mechanisms, A and B, can 
lead to qualitatively new effects. The reason for this is 
very simple: the effectiveness of scattering mechanism A 
in the presence of a current is strongly influenced by the 
momentum dependence of the out-of-equilibrium distri- 
bution of the quasi particles, which is affected by mecha- 
nism B. Therefore, we expect that pab ^ pa + pb even 
in the weak-coupling limit of a semi-classical theory de- 



scribed by a Boltzmann equation (px is the resistivity in 
the presence of scattering mechanism X). 

In this paper we have analyzed such a situation in de- 
tail. Near the AFM quantum-critical point of slightly dis- 
ordered metals, weak isotropic impurity scattering com- 
petes with the strongly anisotropic scattering from spin 
fluctuations, which is most effective along "hot" lines on 
the Fermi surface. The observation that even in rela- 
tively clean single crystals the effective resistivity expo- 
nent a, Ap ot,|T" depends on the amount of disorder 
in the systemQCm, is an important indication that the 
interplay of scattering mechanisms is relevant in these 
systems. Our results for B = explain the experimen- 
tal trends: Ap ex T 15 for dirty systems and Ap ex T 
in high purity samples. A precise test of the spin- 
fluctuation picture which underlies our calculation would 
be to check our scaling predictions using a comparison of 
several high-purity samples with different residual resis- 
tivities po directly at the QCP: Ap « T 1 / 5 /(T/p ) with 
f(y -► 0) = const, and f(y -> oo) ex l/^/y. 

We furthermore propose to use orbital effects in a 
large magnetic field to study anisotropic scattering mech- 
anisms in detail. Interesting nonlinear effects are ex- 
pected if the magnetic field is strong enough to push 
quasi particles within the scattering time Tk over any 
sharp structure with width Afc on the Fermi surface, i.e. 
for e(v x B)rk > Afc. In the limit of very small but 
finite disorder, we were able to calculate the non-linear 
magnetoresistivity analytically in all relevant limits. The 
main effect of a magnetic field is that sharp features in 
the quasiparticle distribution are smeared out. In this re- 
spect both an increase of disorder and a magnetic fields 
lead to similar effects. In a sufficiently strong magnetic 
fields, we therefore expect a T 15 rise of the resistivity at 
the QCP and a T 2 law on the Fermi-liquid side of the 
phase diagram even in samples where the resistivity is 
close to linear in the absence of the magnetic field. Pre- 
cise scaling predictions can be used to check this mecha- 
nism experimentally. 

We hope that the results of this paper establish a set of 
predictions which can be used to investigate mainly one 
question: Can the non- Fermi liquid behavior in clean sys- 
tems be explained by a nearly antiferromagnetic Fermi 
liquid? At the moment this question is still open. Our 
interpretation of the data e.g. of CePd2Si2 favors a SDW 
interpretation. Certainly it is necessary to look also at 
other quantities. For example in CePd2Si2 the Neel tem- 
perature seems to rise linearly with p — p c while in spin- 
fluctuation theory an exponent 2/3 is expected. 

It is important to stress that at least in some systems 
close to an AFM QCP a description within the three- 
dimensional SDW scenario seems to be impossible. The 
best studied example for this is CeCu5.gAuo.1a where 
the spin-fluctuations spectrum shows both an anomalous 
momentumtS and energyta dependence. The linear tem- 
perature dependence of the resistivity in CeCus.gAuo.i 
has a different origin than the one discussed in this paper. 
As the residual resistivity is relatively high, the resistiv- 
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ity should rise with T 15 within the 3d models discussed 
in this paper (a linear resistivity could repilt from 2d 
critical fluctuations coupled to 3d electronsES) . Interest- 
ingly, an exponent 1.5 has been reported for the quantum 
critical resistivity of CeC^-zAg^tS and pressure tuned 
CeCusAuEJ, which might indicate that small changes can 
induce a more conventional SDW-transition in this com- 
pound. CeCu5.9Auo.i1 is Jay, now the most likely and 
best studied candidateHEJo for a trubz_.unconventional 
QCP (YbRh 2 Si2E3 and U .2Yo. 8 Pd 2 Al 3 El have very sim- 
ilar thermodynamic and transport properties). The ce-, 
ported scaling of the susceptibility in CeCu5.9Auo.1t2l 
seems to contradict any simple SDW modelsEa. An in- 
teresting question is whether the rather strong disorder 
is responsible for the anomalous behavior or whether a 
pure interaction effect destroys the Fermi liquid close to 
the QCP. 
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APPENDIX A: APPROXIMATE VARIATIONAL 
SOLUTION OF THE BOLTZMANN EQUATIONS 

In the absence of the magnetic field the Boltz- 
mann equation is equivalent to the following variational 



Ap[Ak] w ar|Afc| 



B 



problemCd'ED which allows to calculate the qualitative be- 
havior of the resistivity in a wry simple and efficient way. 
Following Hlubina and RiceEZI, we define a functional p 
of $k 



dkdk' 



P[*k] = 



iW$k-$ k <) 2 



4e 2 



(/f(v k n)$ k 



min (Al) 



The physical resistivity is given by the minimum of /o[$k] 
regarded as a functional of $k, the scattering matrix F 
is defined in (@) . Instead of solving the Boltzmann equa- 
tion exactly, we can guess the qualitative structure of 
the solution: the strong scattering will equilibrate the 
distribution function in a region of width Afc around the 
hot lines (see Fig. |l|). In the dirty limit, the distribu- 
tion function is of the well known form <E>£ = vrE and 
as a variational ansatz we assume that $k vanishes (or 
more precisely is equal to l/2(4>£ + $ k± o)) in a region 
of width Afc around the hot lines. It is easy to calculate 
the qualitative structure of Ap as a function of the vari- 
ational parameter Afc in the regime where the disorder 
dominates Ap < po : 



(£ + r + (Afc) 2 )( 4 - d )/ 2 



0(i 2 ) (A2) 



The first term is due to the denominator in (Al) and re- 
flects the fact that the disorder favors a smooth distribu- 
tion function. The structure of the second term can easily 
be guessed (or calculated) . It reproduces the well-known 
results e.g. of MoriyaliZl for Afc = 0, who essentially aver- 
ages the scattering rates F k k' over th e fu ll Fermi surface. 
In the dirty limit, the first term in ( A2) dominates, Afc 
vanishes at the minimum and one obtains the well-known 



Ap: 



F 



for Afc < t 



( t + r )(4-d)/2> 

In the ultra-clean limit (i<t 2 ), Afc is of order fc^ and 

Ap « t 2 for Afc « 1. 

For clean systems there is an interesting intermediate 
regime where t + r < (Afc) 2 < 1, where we obtain a 
non-trivial result by minimizing ( |A2J ) 



Ap m (s-^s-a w ith \/t + r < Afc « (t /x) 5 - d < 1. 

This very simple calculation reproduces the results of 
( pij ) qualitatively. 



APPENDIX B: RELAXATION TIME 
APPROXIMATIONS 

Strictly speaking, a relaxation time approximation is 
not valid near a quantum critical point, where the scat- 
tering is dominated by processes involving a fixed mo- 
mentum transfer Q. The time to equilibrate a quasi- 
particle distribution at the point k on the Fermi surface 
will crucially depend on the quasiparticle distribution at 
the point k + Q. For example if Q is perpendicular to 
the electric field E and some reflection or gliding plane 
(and if the Fermi surface is not too complicated), then Q 
will connect points with the same non-equilibrium quasi- 
particle distribution and the singular scattering will not 
relax the non-equilibrium distribution. In such a sit- 
uation vertex corrections cancel the leading self-energy 
contributions and the relaxiation time approximation is 
completely wrong (see also [ |37j on the role of vertex 
corrections for the Hall effects). 

Nevertheless, the relaxation time analysis gives qual- 
itatively the correct resistivity for generic situations, 
where the Q-vectors are not perpendicular to E as can 
be verified by comparison with the full solution of the 
Boltzmann equations discussed in the main text. The re- 
laxation time approximation can have much larger errors 
for the Hall effect because the distribution function <I>k 
perpendicular to the current is often dominated by quasi- 
particles which are scattered into this region by spin- 
fluctuations with momentum Q^ (in the geometry show 
in Fig. |l] this effect is absent) . 
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In second order perturbation, theory the life-time of the 
electron at T = is given byE3 

— = 2^/rV f dwlmx k _ k ,(w)*(w - (e k - e k ,)), (Bl) 
^k j?J 

where x is given by (||) with w q = (q ± Qi) 2 /g 2 . Split- 
ting the integration over k' in a surface integral over the 
Fermi surface and an energy integral perpendicular to it 
J d 3 k' = ff dk'/vp J dew, where Vp(k!) is the Fermi ve- 
locity, performing first the e k ' and then the u) integration 
we obtain 
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^ = r+(5k/ go ) 2 (B3) 



where <5k is the distance of k + Q from the Fermi surface, 
i.e. essentially the distance of k from lines on the Fermi 
surface, where e k = e k ±Q i = p,. Along these "hot lines" 
the inelastic scattering is strongest. 

The scattering rate depends crucially on the distance 
from the hot lines 5k = Sk/kp ~ Sk/qo. Far from the hot 
lines or at lowest energies for a system not directly at the 
quantum critical point, we recover the usual l/r k oc e 2 .. 
Directly at the QCP and for 5k — the scattering rate 
is linear in the quasiparticle energy. At finite tempera- 
tures a typical quasi particle has the energy T and for 
qualitative estimates we use the approximation 
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I I t 2 /(r + t) for 5k < s/r + t 
t 2 /(5 K ) 2 for 8k > s/TTt 



TM 



(B4) 



where tm is a typical scattering rate at the temperature T 
and t = T/T the dimensionless temperature. For r <C 1 
the scattering is highly anisotropic and strongest in a 
region of width 5k rj yr + 1. In the presence of weak 
disorder, the elastic scattering due to impurities has to 
be added 
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where me measure the relative strength of impurity scat- 
tering by the dimensionless quantity x ~ tm/t c \ k, 
1/RRR which is defined in (||). 

The conductivity in relaxation time approximation is 
proportional to the average of r k over the Fermi sur- 
face (assuming that the Fermi velocity does not vary too 
strongly). By integrating (B4) over 5k, we therefore ob- 



tain an estimate for the conductivity. The conductivity a 
in units of the conductivity om = 1/pia is approximately 
given by 
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In the ultra-clean limit, x <JC t < 1, the main contri- 
bution arises from regions far from the hot line or large 
5k. The conductivity diverges with 1/T 2 with a prefactor 
of order 1 



p(Tt) w t 2 for — < min 

T c \ 



1 



or x < t 



(B8) 



as has been pointed out by Hlubina and RiceEZI. In the 

dirty limit, x ^S> t 2 /(r + t), the elastic scattering domi- 
nates and one can expand (B6) or (B7) in the inelastic 
scattering rates: 



Ap(T) 
Pm 



TM 



1 



t z 



for — 3> max 

Tel 



1 



(B9) 
or t <C max[i/rx, x] 



oc 



T 3 / 2 near the 



and we recover the well known 5p 
quantum-critical point. 

In the relaxation time approximation we obtain sim- 
ilar to our results for the Boltzmann equation a pro- 
nounced intermediate regime, where min[l/r^] < 1/t c \ < 
max[l/r^]. This corresponds to the regime where the 
first arctan in (B7) is 7r/2, while the second is small, i.e. 
max[y / rx, x] < t < yfx. The rise of the resistivity in 
linear in temperature in this regime: 

Ap(T) it _ . r /— i r- 

~ t-^v x i° r maxl yrx, x\ <t< \/x. 

Pm 2 



APPENDIX C: NUMERICAL METHODS FOR 
SOLVING THE BOLTZMANN EQUATION 

The Boltzmann equations were solved by direct matrix 
inversion. The Hall effect and the i? 2 -magnetoresistivity 
were calculated in perturbation theory in B. To keep 
the size of the matrices small it is important to use all 
the symmetries of the problem and to choose a suitable 
discretization of the Fermi surface. At low temperatures 
the distribution function $ k changes very rapidly per- 
pendicular to the hot lines. Therefore we define the 
distribution function on polygon-shaped patches on the 
Fermi surface which are relatively long parallel to the 
hot lines Afen ^C kp but short in the perpendicular di- 
rection close to the hot lines Afc <C kpimn[^/i,t/y/x]. 
Because the kernel of the integral equation is strongly 
peaked it is essential to calculate a small number of ma- 



trix elements My = J atch . dk J 



V patch j 



dk'G 



kk> numeri- 



cally if they involve a momentum transfers k — k' close 



1G 



to Qi. To speed up the calculations one can approximate 
I[y] « ir 2 (2n-\-y)/(y(4:Tr 2 +6Try-\-3y 2 )) which is asymptot- 
ically exact in next-to-leading order for large and small 
arguHaejits; errors are smaller than 1.5% (The previously 
usednJ'E3 approximation I[y] rs n 2 /(y(3y + 2n) is also 
sufficient to reproduce next to all of our results with es- 
sentially neglegible error, one exception is the tempera- 
ture dependence of the Hall effect in the ultra-clean limit, 
where it wrongly predicts Rh ~ C\ ± cr{T instead of the 
correct Rh ~ c\ ± C2T 3 / 2 ). We used Fermi surfaces with 
cubic symmetry. As the largest irreducible representation 
of the 48-dimcnsional cubic group is only three dimen- 
sional, it was possible to reduce the size of the matrices 
by the huge factor 16 x 16 using the full symmetry of the 
scattering matrix. We didn't try to solve the Boltzmann 
equations numerically in a finite magnetic field, which 
would break the cubic symmetry, but used perturbation 
theory in B to calculate the figures of this paper. 

It is essential to use a proper discretization of the mag- 
netic field term which conserves 

If [B • (Vk x V k $ k )](rf 2 k/i; k ) (CI) 

I $ k (B • dk) 

JdS 



(d^k x V k $ k ) • B 



for an arbitrary part S of the Fermi surface. Otherwise 
the discretization leads to large errors. To determine 



the averaged v x B contribution for a single patch, (CI) 
is applied for a single patch. $ k on the edges of the 
patch i is determined by a linear interpolation: [B • (v k x 
V k $ k )](d 2 k/ Wk ) « £ j (B.Akj)($ k .,, 2 +$ k .,, 2 )/2 where 
Ak* = k ij]2 — kij-,1 describes the jth edge with corners 
ki,j.2 and kj,j,i- The vectors Ak* point counterclockwise 
around patch i with J] . Ak* = 0. The value of § k at the 
corners ^-i,j, n (n = 1, 2) is determined by a simple average 
over the value of $ on each neighboring patch. The de- 
tails of the se a veraging procedures are not so important 
as long as (|Cl|) is not violated. 



APPENDIX D: MODELS USED FOR THE 
NUMERICAL CALCULATION 

For our numerical calculations we use three differ- 
ent models. Model A is based on a spherical Fermi- 
surface (see Fig. |lj) and the ordering wave-vector is 
Q = (0,0, ±2fcFCOs8# with Qh = 7r/6. The precise 
value of 8h are not very important as long as one stays 
away from 8h = ("2/cf" ordering) or 9h — 7r/2 (fcrra- 
magnetic ordering) where our approach breaks downEj; 
in our numerical calculations we use Qh = 7r/6. The 
dispersion w q in (g) is given by (q — Q) 2 /q 2 . 

Model B describes a cubic system with commensu- 
rate spin waves with ordering vectors Q = (±1,±1,±1) 
and w q = 3 + J2i= x y % cos irqi- The bandstructure was 
choosen to be featureless but it includes terms which 



mix the various directions and it is characterized by 
8 non-intersecting "hot lines" : e k = t J2 X y z cos n h + 
t' J2 X v z cos 27r/ci + i"(cos Trk x cos Trk y + cos Trk x cos irk z + 
cos nky cos nk z ) - \i with t = -0.2, t' = 0.05, t" = -0.045 
and fi = -0.0225 

We have compared model B to model C, where the 
ordering vector is parallel to the principal axes: Q; = 
(±1, 0, 0), (0, ±1, 0), (0, 0, ±1) using also a different band- 
structure. All results are very similar compared to those 
of model B and are therefore not shown in this paper. 
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